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Abstract 


Nowadays, the creation of a spatial distribution map for assessing fire risk is widely recognized 
as a crucial management tool at various levels. This tool helps monitoring natural resource 
sustainability and the effective control of this environmental hazard. Through the integration of 
field operations, remote sensing data, geographic information system techniques, and diverse 
statistical methods, it becomes feasible to develop a dependable spatial prediction of fire hazard 
potential for different regions. In this research, nine factors were identified to be effective in fire 
risk modeling, including altitude, slope, aspect, distance from the road, NDVI, LST, TWI, and 
TPI. Fuzzy Analytic Hierarchy Process (FAHP) and Logistic Regression (LR) were employed to 
identify risk areas and determine the most significant factors influencing the occurrence and 
spread of fire. Historical fire areas were identified using Google Earth Engine and MODIS 
images. The initial results showed that both models assign the highest coefficients to NDVI, LST, 
and distance from the road. However, during the verification phase, the performance characteristic 
curve of both models was relatively similar (0.847 for FAHP method and 0.837 for LR method). 
Upon examining historical fire pixels, it was found that FAHP method correctly identified 
approximately 87% of the pixels belonging to classes with high and very high risk, whereas LR 
method only overlapped with 2246 of these pixels. This suggests that FAHP method is better at 
identifying areas with a high risk potential compared to LR method. While it is important to 
acknowledge that the creation of risk prediction maps using various models cannot completely 
eliminate all fires, it can greatly diminish their frequency and help their management by offering 
effective solutions. 
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مقاله پژوهشی 
ارزیابی پتانسیل خطر آتش‌سوزی با استفاده از رویکردهای تحلیل سلسله مراتبی فازی و رگرسیون لحستیک 


0 ہی تیموری" - استادیار گروه مهندسی طبیعت. دانشکده کشاورزی شیروان. دانشگاه بجنورد. بجنورد. ایران. 
آیدینگ کرنژادی- استادیار پژوهشی» بخش تحقیقات منابع طبیعی مرکز تحقیقات و آموزش کشاورزی و منابع طبیعی استان گلستان» سازمان 


تحقیقات. آموزش و ترویج کشاورزی, گرگان ایران. 


تاریخ دریافت: ۱۶۰۱/۹/۲۰ تاریخ بازنگری: ۱۶۰۱/۱۱/۲۵ تاریخ تصویب: ۱۶۰۱/۱۲/۸ 


حکده 


امروزه تهیه نقشه توزیع مکانی خطر آتش‌سوزی یکی از ابزارهای ضروری مدیریت در سطوح مختلف 
جهت پایش پایداری منابع طبیعی و JES‏ این مخاطره محیط‌زیستی است. تلفیق عملیات میدانیء داده‌های 
دورسنجی, تکنیک‌های سیستم اطلاعات جغرافیایی و آماری مختلف می‌تواند پیش بینی فضایی قابل‌اعتمادی 
از پتانسیل خطر آتش‌سوزی برای مناطق مختلف ایجاد کند. در این تحقیق از ٩‏ عامل مؤثر در مدلسازی 
خطر آتش‌سوزی شامل عوامل ارتفاع شیب جهت. انحنای دامنه. فاصله از جاده و شاخص‌های NDVI‏ 
TWI LST‏ و TPI‏ با دو روش تحلیل فازی سلسله مراتبی و رگرسیون لجستیک برای شناسایی مناطق 
خطر و تعیین مهم‌ترین عوامل مؤثر در بروز و گسترش آتش استفاده شد. برای شناسایی مناطق حریق 
تاریخی نیز از سامانه ابری گوگل ارث انجین و تصاویر مادیس استفاده شد. نتایج اوليه OUS‏ داد در هر دو 
مدل» عوامل LST NDVI‏ و فاصله از جاده بیشسترین ضرایب را به خود اختصاص داده‌اند. در مرحله 


صحت سنجی نیز اگرچه منحنی مشخصه عملکرد هر دو مدل» نسبتاً یکسان (مقدار ۰/۸۶۷ در روش تحلیل 
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سلسله مراتبی فازی و مقدار ۰/۸۳۷ در روش رگرسیون لجستیک) بود. با بررسی پیکسل‌های حریق 
a6 ja‏ ساس د ور وق تال سال ی eR‏ دوو درد گا هاش قات شظر 
زیاد و خیلی 285 و در روش رگرسیون لجستیک فقط ۲۲ درصد پیکسل‌های طبقات مزبور با مناطق دارای 
ساق ریق Tk‏ اق للا کر مرش رو فطل مقا E‏ ی توف 
رگرسیون لجستیک توانسته مناطق با پتانسیل خطر بالا را شناسایی کند. اگرچه تهیه نقشه‌های پیش‌بینی 
خطر توسط مدل‌های مختلف از وقوع کلیه حریق‌ها جلوگیری نخواهد کرد ولی می توان با ارائه راهکارهای 
مدیریتی, وقوع آن را کاهش داده و کنترل OT‏ را تسهیل نمود. 

کلیدواژه‌ها: تصمیم گیری چندمعیاره» پیش بینی خطر. گوگل ارث انجینء منحنی مشخصه عملکرد. ماهواره 


مادیس. 


۱- مقدمه 


جنگل‌ها یکی از مهم‌ترین منابع طبیعی تجدید شونده بوده و به‌عنوان سوپاپ اطمینان» تعادل اکولوژیکی جهانی و 
پایداری تمدن بشر را حفظ می کنند (نوتاماکت و استراتولیاس ؛ ۲۰۲۱). اگرچه معیارهای مختلفی با توجه به نوع اقلیم 
و شرایط جغرافیایی منطقه به‌عنوان شاخص‌های مدیریت پایدار منابع طبیعی وجود دارد. لیکن اصول اولیه تمام آن‌ها 
پر اساس سلامت اکوسیستم» تجدید حیات عدم روند تخریب و عدم‌تغییر کاربری آن استوار است (زارع چاهوکی؛ 
۸ یکی از مهم‌ترین مصادیق تخریب و تهدید Ol‏ آتش‌سوزی سیستم‌های جنگلی در سطح جهان است 
(جعفری و همکاران» ۲۰۱۹). وقایع آتش‌سوزی اثرات منفی بر تنوع زیسستی بیلان اکولوژیکی. حیات‌وحش؛ 
فرایندهای فیزیکی و شیمیایی زمین. تغییر اقلیم» خصوصیات هیدرولوژیکی فعالیت‌های انسانی و اجتماعی-اقتصادی 
دارد (ونکاتش" و همکاران» ۲۰۲۰). اگرچه مهار کامل این پدیده امکان‌پذیر نیسست ولی با اقداماتی می‌توان به 
پیشگیری بخشی و کاهش اثرات سوء آن دست یافت. یکی از اقدامات ممکن, پهنه‌بندی خطر از نظر احتمال وقوع و 
تعیین مناطق با پتانسیل بالای حریق است (اسکندری و اسکندری. Mert‏ به‌صورت کلی» چهار عامل توپوگرافی, 
اقلیمی» پوشش گیاهی و فعالیت‌های انسانی به‌عنوان عوامل اشتعال‌زا در جنگل‌ها محسوب می‌شوند لاک و 
سونار" ۲۰۲۰). همچنین روش‌های مختلف مدل‌سازی خطر آتش‌سوزی شامل رگرسیون خطی, رگرسیون لجستیک» 
تحلیل تصمیم چند معیاره و روش‌های همپوشانی وزنی مورداستفاده قرار گرفته است (عابدی. ۲۰۲۲؛ فرامرزی و 
همکاران» ۲۰۲۱). روش فرایند تحلیل سلسله مراتبی به‌عنوان یکی از پرکاربردترین روش‌های مکانی با موفقیت بالا 
در شناسایی چنین مناطقی در بسیاری از مناطق دنیا مورداستفاده قرار گرفته است (پورقاسمی و همکاران .)۲۰۱٦‏ مثلا 


1 Nuthammachot and Stratoulias 
2 Venkatesh 
3 Colak and Sunar 
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در ایرانء عابدی قشلاقی و همکاران (۱۳۹۷) در تحقیق خود به ارزیابی و پھنەبندی خطر آتش‌سوزی جنگل با استفاده 
از تکنیک‌های تصمیم‌گیری چند معیاره و ساج در جنگل‌های نوشهر پرداخته و نتیجه گرفتند که حدود ۲۲ درصد 
منطقه در کلاس خطر 265 و بسیار زیاد قرار دارد. تومار" و همکاران (۲۰۲۱) نیز به ارزیابی خطر حریق در جنگل‌های 
هند با استفاده از تحلیل سلسله مراتبی و مدل‌سازی مکانی پرداخته و مناطق با طبقات مختلف حریق را مشخص 
نمودند. بوسیکو" و همکاران (۲۰۱۹) نیز به مدل‌سازی آتش تحت شرایط تغییر اقلیم در جنگل‌های ایتالیا پرداخته و 
نتایج آن را قابل‌قبول دانستند. 
در سال‌های اخیر تلفیق روش‌های تصمیم‌گیری نوین چندمعیاره در گستره مکانی و همچنین مدل‌های آماری 
جدیدت چارچوب مناسبی را به‌منظور ارزیابی مناطق مستعد آتش‌سوزی قرار داده است (نوو" و همکاران, ۲۰۲۰). 
مثلاً اردین و چاگلار* (۲۰۲۱) و عابدی قشلاقی و همکاران (۲۰۲۰) به ارزیابی خطر آتش‌سوزی با ترکیب روش 
تحلیل سلسله مراتبی با منطق فازی پرداخته و نتایج حاصله را منطقی و قابل‌قبول‌تر نسبت به روش تحلیل سلسله 
مراتبی دانستند. یاتیش" و همکاران (۲۰۱۹) به استفاده از دو روش رگرسیون لجستیک و روی‌هم‌گذاری وزنی در 
ترسیم مناطق خطر جنگل‌های هند پرداخته و روش رگرسیون لجستیک را روش مناسب‌تری بخصوص در شرایط 
وجود داده تاریخی عنوان نمودند. 
یکی از ملزومات مدل‌سازی آتش» دسترسی به داده‌های باکیفیت زمانی و مکانی مناسب است (پیائو و همکاران, 
۲ ماهواره مادیس این قابلیت را همراه با ارائه محصولات متنوع پردازشگر فراهم آورده است. مثلاً شریف‌نژاد و 
همکاران (۱۶۰۰) به ارزیابی قابلیت محصولات آتش سنجنده مادیس در استان گلستان پرداخته و اذعان نمودند که 
این محصولات فقط نیمی از حریق‌های صورت گرفته را OUS‏ دادند و در مورد آتش‌سوزی‌های جزئی ضعیف عمل 
می‌کنند که این موضوع را به مقیاس مورداستفاده این محصولات که در سطوح بزرگ کارایی دارد مرتبط دانستند. 
هیسلوپ" و همکاران (THY)‏ نیز به بررسی سری زمانی داده‌های حریق جنگل‌ها در تصاویر مادیس و شاخص 
نسبت نرمال آتش‌سوزی پرداخته و نتیجه گرفتند که بین سال‌های ۲۰۰۱ تا ۲۰۱۸ء حدود ۱۸۰ میلیون هکتار از Y‏ 
میلیارد هکتار جنگل‌های مناطق معتدل مورد مطالعه» دچار حریق شده و روند افزایش حریق در جنگل‌های شیلی و 
روند کاهشی در جنگل‌های چین مشاهده می‌شود. 
هرچند تصاویر ماهواره‌ای مثل مادیس کمک شلیانی در این امر می‌کنند. حجم بالا و زمان طولانی پردازش 
محدودیت‌هایی را به همراه دارد. لذا استفاده از محصولات ماهواره‌ای که دارای قابلیت تحلیل بدون ذخیرہ داده‌ها 
Tomar‏ 1 
Busico‏ 2 
Novo‏ 3 
Erdin and Çağlar‏ 4 
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باشند بسیار مفید به نظر می رسد (پیائو و همکاران. ۲۰۲۲). این موضوع توسط سامانه پردازش ابری گوگل ارث 
انجین مهیا شده است. چاوداراوغلو! (۲۰۲۱) و glo‏ و همکاران (۲۰۲۲) به مدت‌زمان پردازش سریع این سامانه در 
تشخیص مناطق Gao‏ اشاره داشته‌اند. 

بررسی سوابق مطالعاتی گذشته در داخل و خارج از olal‏ نشان می‌دهد که جهت تعیین پتانسیل خطرات 
آتش‌سوزی‌هاء استفاده از داده‌های میدانی و مدل‌سازی مکانی توأم با داد‌های تاریخی ماهوارهای که قابلیت تحلیل 
سری زمانی حریق را فراهم می‌سازد بسیار مناسب هستند. لذا در این تحقیق از قابلیت مدل‌سازی فازی در تحلیل 
سلسله مراتبی و مقایسه آن با مدل رگرسیون لجستیک همراه با صحت‌سنجی تصاویر حریق ماهواره مادیس در محیط 
گوگل ارث انجین برای تعیین مهم‌ترین عوامل مؤثر در گسترش و تهیه مدل و نقشه خطر آتش‌سوزی در پارک ملی 
گلستان استفادہ شدہ است. 


-Y‏ مواد و روش‌ها 

۱-۲- منطقه مورد مطالعه 

پارک ملی گلستان به‌عنوان قدیمی‌ترین پارک ملی ثبت شده در ایران در محدوده Y‏ استان گلستان و خراسان 
شمالی واقع است. محدوده پارک بین طول‌های جغرافیایی ٩۳‏ 1۳ 00 تا EY"‏ 9۵۱۱۷ شرقی و عرض‌های 
جغرافیایی ۲۱۳ vv‏ ۰۳۷ تا YA oa"‏ 0۳۷ شمالی eil‏ است و مساحتی حدود AOA‏ کیلومترمربع و محیط ۱۶۵ 
کیلومتر دارد. اقلیم منطقه در بخش‌های مرتفع» معتدل کوهستانی و در ارتفاعات کمتر اقلیم معتدل و نیمه مرطوب 
است. در شکل ۱ موقعیت پارک ملی گلستان در ایران و استان‌های خراسان شمالی و گلستان نشان داده شده است. 


1 Cavdaroglu 
2 Piao 
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کت سی 


شکل 4- موقعیت پارک ملی گلستان در ایران و استان‌های خراسان شمالی و گلستان 


-Y-Y‏ داده‌های مورداستفاده 

ازآنجایی که ماهواره مادیس, اولین ماهواره‌ای است که از سنجنده‌های حرارتی برای پایش آتش‌سوزی‌ها و ثبت 
تصاویر آن‌ها استفاده می‌کند KE)‏ لیا" و همکاران ۲۰۰۳) لذا از این ماهواره برای شناسایی نقاط آتش‌سوزی ضمن 
تطبیق با گزارش‌های اداره کل منابع طبیعی و آبخیزداری استان استفاده شد. برای پردازش داده‌های مزبور از سامانه 
پردازشی ابری گوگل ارث انجین استفاده شد. شاخص مورداستفاده در شناسایی آتش‌سوزی‌ها شاخحص نسبت 
سوختگی نرمال 0 است (معادله 6 

NBR = (NIR — SWIR)/(NIR + SWIR) (۱) 


که NIR‏ باند مادون‌قرمز نزدیک و SWIR‏ باند مادون‌قرمز طول موج ab S‏ است. 


1 Giglioa 
2 Normalized Burn Ratio (NBR) 
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با گذشت مدتزمانی از آتش‌سوزی و رشد مجدد/ بازسازی پوشش گیاهی» شاخص NBR‏ تأثیر کمتری دارد. لذا 
شاخحص dNBR‏ به‌عنوان ابزاری برای درک بهتر وسعت و شدت آتش, پس از محاسبه تفاوت بین شرایط قبل و بعد 
از آتش‌سوزی (معادله ۲) ارائه شد (فاسناچت ' و همکاران ۲۰۲۱). 
dNBR = PrefireNBR - (Y)‏ 
با توجه به لحاظ نمودن موقعیت مکانی. زمانی» مطالعات قبلی و نظرات کارشناسان محلی, تعداد ۹ عامل شامل 
ارتفا شیب جهت. فاصله از جاده. انحنای دامنه و همچنین شاخص‌های TWP 1۳1 010۷1 LST?‏ که به ترتیب 
Sl‏ دمای سسطحی, تراکم پوشش گیاهی. موقعیت توپوگرافی و رطوبت منطقه هستند به‌عنوان عوامل مؤثر در 
شروع و گسترش آتش‌سوزی جهت استفاده در مدل‌های تحقیق انتخاب شد. 
ازآنجایی که نوع پوشش گیاهی در هر منطقه معمولاً بەصورت ترکیبی از جوامع گیاهی مختلف است. لذا استفاده 
از شاخص‌ھای پوشش گیاهی مستخرج از تصاویر ماهواره‌ای بسیار معمول و کاربردی است. یکی از این شاخص‌هاء 
شاخص NDVI‏ است که بیانگر فرایند فتوستز و درنتیجه محتوای رطوبت گیاه و خاک است (لی" و همکاران, 
YY‏ که در تصاویر لندست ۸ از معادله ۳ قابل‌تعیین است. 
NDVI = (NIR — R)/(NIR + R) (۳‏ 
که در رابطه فوق» UL NIR‏ مادون‌قرمز نزدیک (۰/۷۲-۰/۹ میکرومتر) و R‏ محدوده باند قرمز (۰/۱۳-۰/۷۹ 
میکرومتر) است (یاتیش و همکاران ۲۰۱۹). 
همچنین ازآنجایی که پتانسیل حریق با کاهش رطوبت دامنه افزایش پیدا می‌کند. استفاده از شاخصی که بیانگر 
رطوبت مسطحی و کنترل کننده گسترش آتش اسست بسیار حائز اهمیت اسست (ماسیندا" و همکاران ۲۰۲۱ 
معمول‌ترین شاخصی که به تمایل پیکسل‌ها در نگهداشت رطوبت می‌پردازد شاخص TWH‏ است (ماتیوی و همکاران, 
۵۹ لذا در این تحقیق از شاخص مزبور (معادله )٤‏ استفاده شده است. 
TWI = In(a/tanB) (£)‏ 
که » مساحت تجمعی زھکشی بالادست و tanp‏ زاویه شیب هر نقطه است که در این مطالعه برای تهیه فایل 
آن از نرم‌افزار 5464-615 (کنراد" و همکاران, ۲۰۱۵) استفاده شده است. مناطق مستعد به تجمع آب (شبکه زھکشی) 
CYL polis‏ شاخحص و مناطقی که زهکشی خوبی دارند مقادیر کمتری دارند. اگرچه بررسی دقیق این شاخص به 
Fassnacht‏ 1 
Land Surface Temperature‏ 2 
normalized difference vegetation index‏ 3 
Topographic Position Index‏ 4 
topographic wetness index‏ 5 
Li‏ 6 


7 Masinda 
8 Conrad 
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نوع اقلیم» نوع دامنه و مقیاس Of angi‏ بستگی دارد (آگرن" و همکاران, ۵۶ از سایر متغیرهای اقلیمی مهم در شیوع 
و گسترش آتش نیز می‌توان به دمای سطح زمین اشاره کرد که برای بررسی دقیق آن از شاخص LST‏ و محصول 
آمادہ MOD11A2‏ که داده‌های حداقل و حداکثر دمایی را بەخوبی ثبت می کند (آیانلادہ', )۲۰۱٦‏ استفاده می‌شود. برای 
بررسی داده‌ها و شاخص‌های Sd‏ شدہہ تمام عوامل در قالب فایل رستری با اندازه پیکسل ۳۰ متری تهیه شدند. 

-Y-Y‏ روش AHP‏ فازی 

۳ یک رویکرد تصمیم گیری چند معیارہ است که در آن معیارها در یک ساختار سلسله مراتبی سازمان‌دهی 
می‌شوند (زارعیان و همکاران. ۲۰۲۱). در این رویکرد اهمیت نسبی معیارها بر اساس وزن آن‌ها L)‏ توجه به نظرات 
کارشناسان» متخصصان و بررسی سوابق قبلی) مشخص شده و ماتریس مقایسات زوجی ساخته می‌شود که 
نشان‌دهنده میزان قوت هر معیار (از عدد یک تا )٩‏ نسبت به سایر معیارها است. سپس این مقادیر نسبی, به کمیت‌های 
عددی تبدیل شده و درنهایت باید بزرگ‌ترین مقدار ویژه Amar)‏ در محاسبه شاخص ناسازگاری (CI)‏ محاسبه شود 
(معادله 0( 

CI = (Amax — n)/(n — 1) (6) 

که 7 تعداد معیارهای مورد بررسی است. نرخ ناسازگاری نیز با مقایسه شاخص ناسازگاری با شاخص تصادفی 
که توسط Y gela‏ (۱۹۸۶) ارائه شده از معادله ٦‏ قابل محاسبه است. 


CR = CI/RI C) 
چنانچه عدد مزبور کوچکتر یا مساوی ۰/۱ باشد سازگاری سیستم قابل‌قبول است وگرنه باید در قضاوت‌ها‎ 
تجدیدنظر نمود.‎ 


گاهی col‏ هنگام ایجاد یک ماتریس مقایسه نمی‌توان مقادیر واضحی را برای مقایسه معیارهای مختلف تین 
کرد. در این موارد می‌توان از منطق فازی بهره گرفت که اعتبار هر معیاری را با طیف خاکستری (عدد بین صفر و 
یک) نشان می‌دهد. لذا با اعمال تابع عضویت خطی و جمع لایه‌های وزن‌دان نقشه نهایی خطر آتش‌سوزی به‌صورت 
رستری تهیه می‌شود. بدیهی است هرچه مقدار پیکسلی به عدد یک نزدیک شود Sly‏ خطر بالاتر است. همچنین 
نقشه تهیه شسدہ با طبقه‌بندی مجدد و به روش جنکز (آجین * و همکاران» ۲۰۱۷) قابل‌تقسیم به طبقات مختلف 


1 Agren 

2 Ayanlade 
3 Saaty 

4 Ajin 


VA‏ جغرافیا و مخاطرات محیطی شمارۂ سوم 


۶-۲-روش رگرسیون لجستیک 


رگرسیون لجستیک یک روش از مدل‌های موسوم به خطی تعمیمیافته است که قابلیت پیش بینی خروجی‌های 
مجزا از مجموعه‌ای an gy,‏ یا غیریبوسته را دارد. استفاده از این روش در مجموخه داده‌های غیرعطی بسیار کارا برده 
و هدف اصلی آن یافتن بهترین مدل برای بیان رابطه بین یک متغیر وابسته و چندین متغیر مستقل است (چویکو' و 
همکاران. ۲۰۱۰). معادله اساسی مدل مزبور به‌صورت معادله V‏ است. 

In(p/(1 — p)) = Bo + بترم +۰۰۰4 + درم‎ (V) 

که D‏ احتمال وقوع پیکسل مکانی آتش‌سوزی, X‏ عوامل مؤثر و B‏ ضرایب اکتشافی متغیرهای معادله هستند 
(حیدری و همکاران» ۱۶۰۰). در تحقیق حاضس متغیر وابسته وقوع یا عدم وقوع آتش‌سوزی و عوامل ٩‏ گلنه 
متغیرهای مستقل هستند. 


ROC منحنی‎ —-0o-Y 


ازانجایی که اعتباریابی یک فرایند ضروری در ارزیابی مدل‌های پیش‌بینی است. از روش منحنی مشخصه عملکرد 
سیستم (ROC-AUC)‏ که یکی از متداول‌ترین روش‌های کمی برای پهنه‌بندی و ارزیابی مخاطرات طبیعی امستفاده 
می‌شود که انعطاف‌پذیری زیادی در بیان درجه حساسیت مدل‌ها دارد eld)‏ و همکاران» ۲۰۲۱). این منحنی بر اساس 
تقسیم درصد مثبت‌های صحیح بر درصد منفی‌های کاذب ترسیم می شود و سطح زیر آن بیانگر احتمال پیش‌بینی 
صحیح یک پیکسل تصادفی است (ابراهیمی و همکاران, ۱۳۹۷) و به‌صورت محدوده‌های ۰/0-۰/۷ AIAN‏ ۰/۸- 
۷ ۰-۸۹/ و ۰/٩-۱‏ به ترتیب پیش بینی ضصعیف. متوسط خوب. خیلی خوب و عالی دسته‌بندی می‌شود 
(پورقاسمی. ۲۰۱۵). 
۳- نتایج 

عوامل ٩‏ گانه موثر بر آتش‌سوزی که در این تحقیق مورداستفاده قرار گرفتند به‌صورت لایه‌های رستری و با اندازه 
پیکسل‌های یکسان تهیه شسد. سری Shy‏ قاط حریق ہی سال‌های ۲۰۰۶ 1۶ ۲۰۲۰ نیز از ماهوازه مادیس hag‏ 
کدنویسی در سامانه گوگل ارث انجین استخراج (شکل ۲) و پس از انطباق با داده‌های مشاهدات زمینی گستره آن‌ها 
برای اهداف واسنجی و صحت‌سنجی مشخص شد. تعداد ۱۲۲ موقعیت آتش‌سوزی بین سال‌های ۲۰۰۶ تا ۲۰۲۰ 


استخراج شد که از ۰ درصد آن‌ها برای واسنجی و ۰ درصد برای صحت‌سنجی استفاده شد. لایه‌های موضوعی 
Ag‏ شده در شکل Y‏ نشان داده شده است. 


] Chuvieco 
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شکل -Y‏ نقشه مناطق دچار آتش‌سوزی در سال‌های ۲۰۰6 تا ۲۰۲۰ مستخرج از سامانه گوگل ارث انجین در 
پارک ملی گلستان 
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شکل ۳- نقشه‌های موضوعی پارک ملی گلستان شامل نقشه ارتفاعی (a)‏ نقشه شیب (b)‏ نقشه جهت شیب )€( 


(h) TPI نقشه‎ (g) Wetness نقشه‎ (f) LST نقشه‎ de) نقشه انحنای دامنه‎ (d) NDVI نقشه‎ 


همان‌طور که در شکل ۳ ملاحظه می‌شود. ارتفاع پارک ملی گلستان بین EW‏ تا ۲۳۵۹ متر از سطح bya‏ با میانگین 
۵ متر و حداکثر شیب طولی منطقه ۱۹/٤‏ درجه است. جهت دامنه‌ها عمدتاً شمالی (حدود ۳۶ درصد منطقه) و 
قسمت متراکم پوشش گیاهی در بخش غربی پارک و در محدوده استان گلستان قرار دارد. ازآنجایی‌که شاخص 
1 معادل عدد یک بیانگر پوشش WE‏ متراکم است» وجود مقادیر بالاتر از ۰/۵ در بخش بزرگی از منطقه بیانگر 
پوشش جنگلی است. بررسی نقشه‌های کاربری اراضی رقومی شده نیز نشان می‌دهد حدود 1۷/۵ درصد از پارک 
پوشیده از جنگل متراکم» ۳۳ درصد جنگل تنک و مابقی به بیشه‌زارها و مراتع اختصاص دارد. در مورد شاخص 
انحنای دامنه نیز با توجه به اينکه مقادیر مثبت و منفی به ترتیب بیانگر شیب‌های محدب و مقعر اسست. علیرغم دامنه 
تقریبی -٩‏ تا ۹+ این شاحص. مقدار میانگین ۰/۰۲۸- و بررسی هیستوگرام شاخص در هر پیکسل نشان‌دهنده OF‏ 
است که بیشتر دامنه‌ها دارای شیب یکنواختی هستند. بررسی نقشه LST‏ نیز OUS‏ می‌دهد این متغیر تقریباً بین ۲۲/۵ تا 
۷ درجه با میانگین ۳۱/۶ در قسمت‌های مختلف منطقه نوسان داشته و مقادیر کم در قسسمت غربی Sok‏ 


نشان‌دهنده اثر تعدیلی جنگل‌های متراکم بر دما است و طبیعی است که بخش شرقی منطقه که جنگل کم تراکم و 
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پیشه‌زارها را به خود اختصاص می‌دهد دمای بالاتری را نشان دهد. نقشه TWI‏ نیز حداقل مقدار ۲/۵۳ و حداکٹر ۱٦/۷‏ 
را نشان می‌دهد که مقدار عددی متوسط ۷۱۷ بیانگر شرایط زهکشی متوسط است. لایه موضوعی 1۳1 (شکل h‏ نیز 
شاحص موقعیت تویوگرافی را نشان می‌دهد که بین ۱۸۵- تا ۱۸۵+ با میانگین ۱/۰۵- متغیر بوده که نشان‌دهنده این 
است که به‌صورت معمول منطقه دارای شیب‌های میانی. یال‌های موضعی و دره‌های باریک تا شکل است. با توجه به 
گزارش‌های فنی قبلی ازجمله گزارش منتشر شده توسط مرکز پژوهش‌های مجلس (۱۳۹۹) که عامل انسانی را یکی 
از دلایل اصلی بروز آتش‌سوزی‌ها عنوان نموده است. نقشه موضوعی فاصله از جاده و حریق‌های ثبت شده 
به‌صورت نقشه جداگانه در شکل ٤‏ آورده شده است. همان‌طور که ملاحظه می‌شود ST‏ آتش‌سوزی‌های رخ داده در 
حاشیه یا فاصله نزدیکی از جاده‌ها به وقع پیوسته‌اند که احتمال Lice‏ انسانی (عمدی یا غیرعمدی) آن‌ها را تقویت 
می‌کند. بخصوص قسمت جنوبی پارک که جاده آسیایی گرگان- مشهد از بین آن می‌گذرد و میزان توقف و اسکان 
مسافران و گردشگران بسیار زیاد است. 


37*30'0"N 


37°20'0"N 


55°50'0"E ع5600‎ 56°10'0"E 


شکل -٤‏ نقشه پراکنش جاده در سطح Syl‏ و مناطق دچار آتش‌سوزی در سال‌های ۲۰۰6 تا ۲۰۲۰ در پارک 
ملی گلستان 


\AY‏ جغرافیا و مخاطرات محبطی شمارۂ سوم 


پراکنش زمانی آتش‌سوزی‌های مورد بررسی در ماه‌های مختلف در شکل ٥‏ نشان داده شده است. ملاحظه 
می‌شود که در clack‏ تیر الی OUT‏ که مصادف با دمای Mb‏ پوشش گیاهی و زیر اشکوب جنگلی متراکم‌تر و تردد 
بالای گردشگران است. فراوانی آتش‌سوزی‌ها بیشتر است. این موضوع فرض تلفیق دلایل طبیعی و غیرطبیعی 
به‌عنوان Lie‏ بروز و گسترش حریق را تقویت می‌کند. علی نیا و همکاران OE)‏ نیز فراوانی وقوع حریق‌های ماه 
خرداد تا شهریور را در مراتع» جنگل‌ها و اراضی زراعی لرستان گزارش نموده بودند. بروز آتش‌سوزی در ماه مهر و 
آبان نیز با کاهش رطوبت نسبی هوا و وزش بادهای گرم موسسمی مرتبط است که در گزارش مرکز پژوهش‌های 
مجلس (۱۳۹۹) اشاره شده است. 
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شکل ۵- تعداد وقایع حریق رخ داده مورد بررسی در ماه‌های مختلف بین سال‌های ۲۰۰۶ تا ۲۰۲۰ 


برای درک ارتباط بین متغیرهای مورد بررسی و آتش‌سوزی‌های رخ دادہ با برش فایل رستری لایه‌های موضوعی 
با مرز حریق‌های تاریخی در سامانه گوگل ارث انجین» به بررسی متغیرها در پیکسل‌های دچار حریق پرداخته شده 
است که نتایج آن در جدول ۱ ملاحظه می‌شود. برای تبیین تغییر میانگین متغیرها در مناطق دچار حریق نسبت به 
میانگین شرایط پارک از روش شکستگی‌های طبیعی (جنکز) استفاده شده است. بدین ترتیب که دامنه حداقل و 
حداکثر متغیر به ۵ طبقه تقسیم شده و مقایسه صورت می‌پذیرد. همان‌طور که در جدول ۱ مشاهده می شود فقط در 
مورد متغیر ارتفاع و شیب تغیبر طبقه صورت گرفته است و بدین معناست که فقط میانگین دو پارامتر ارتفاع و شیب 
obla‏ فان سرع مرو اگ obo a Logs cnin ot Se‏ سی خاضی فانسائل 

اله ور عائل عم pe ya‏ 6۷1 فال کین با روش کر تس را پررسی لاه یکل cei‏ 
جهت‌های مختلف استنباط شد که کماکان جهت دامنه‌های شسمالی جهت غالب بوده و برای شاخص NDVI‏ هم 


مقادیر میانگین حدود ۰/۸ محدوده پوشش جنگلی و بیشه‌زارهای متراکم است که تغیبری نیافته است. 


سال دوازدهم ارزیابی پتانسیل خطر آتش‌سوزی با استفاده از .... ۱۸۳ 


جدول ۱- مقادیر میانگین. حداقل و حداکثر متغیرهای مورد بررسی در سطح پارک و در محدوده حریق‌های 


تاریخی 
m 7‏ سطح پارک ملی گلستان " سطح حریق‌های مشاهداتی 7 
میانگین حداقل حداکٹر میانگین حداقل حداکثر 
ارتفاع متر ۱۲۵ t YYAs ۸۷۳ 1110 Y ۲۳۹۹ ew‏ 
شیب درجه ۷۷۸۵ VEE‏ ۲ ۷۷۳۱ . 10/7 ۳ 
شمالی (۲۳/۷درصد)- جنوبی (۲۰درصد)- شمالی (۲۸درصد)- جنوبی (۲۲درصد)- 

j غربی (۱۸۷/۶درصد)- شرقی (۱۷/۸درصد) : غربی (۱۹درصد)- شرقی (۱۷درصد)‎ ۲ a 
Y ۵۱ - ۲ qi Y +4 -۹ -۸ z انحنای دامنه‎ 
: ۰/۹۲ ۰۳ /۸ z 4 ۰٦ /۵ z NDVI 
Y M -۹٦ HA Y +۵ ET -۵ z TPI 
Y ٦ ۲/٦٢ vi Y vw yoy KAM 3 TWI 
Y EUN ۷ ۳/۹ ۳ 11۷ ۳۳/۵ ۳/۶ سلسیوس‎ 18 


در روش AHP‏ فازی, بعد از مرور عوامل مؤثر در آتش‌سوزی‌ها و تعیین وزن هر متغیر توسط ale‏ و مستندات» 
ماتریس مقایسات زوجی به‌عنوان ورودی به ابزار تشکیل و نرخ سازگاری با عدد ۰/۰۸ به دست آمد که با توجه به 
اینکه کمتر از ade‏ ۰/۱ است ماتریس به دست آمده قابل‌قبول می‌باشد. وزن‌های به دست آمده از ماتریس مقایسات 
زوجی برای ag‏ معادله خطر آتش‌سسوزی (A ablas)‏ در نرمافزار 100 مورداستفاده قرار گرفت. پیکسل‌های 
خروجی به‌صورت طبقه‌بندی شده در ۵ طبقه استخراج که مساحت هر طبقه در جدول Y‏ آورده شده است. ملاحظه 
می‌شود که حدود ۲۲ درصد پارک در این روش پتانسیل زیادی در بروز آتش‌سوزی داشته و تقریاً نیمی از پارک 


خطر آتش‌سوزی متوسط به پایین EWA)‏ درصد) و نیمی دیگر (۵۲/۲ درصد) دارای خطر متوسط به بالا هستند. 


FRrayp = 0.29670 + 0.18LST + 0.1746N + 0.10295 — 0.0862TW + 0.08234 + 0.0332H + 0.02040 — 
0.0237TP (A) 


FR aS‏ خطر آتش‌سوزی»ظ فاصله از حاده LST‏ دمای سطحی. N‏ شاخص S «NDVI‏ شیب دامنف TW‏ شاخص 
همان‌طور که در معادله A‏ مشخص است. عوامل فاصله از جاده L)‏ ضریب ۰/۲۹)» دمای سطحی (۰/۱۸) و تراکم 
پوشش گیاهی (۰/۱۷) از مؤثرترین عوامل بروز و گسترش حریق در پارک محسوب می‌شوند و ضرایب بیشتری را 


به خود اختصاص داده‌اند. 


ء۸ جغرافیا و مخاطرات محبطی شمارۂ سوم 


جدول ۲- مساحت طبقات خطر آتش‌سوزی در روش فازی سلسله مراتبی 


درجه خطر آتش‌سوزی مساحت (کیلومترمربع) درصد درصد تجمعی 
خیلی VW VW oV/YA e‏ 
کم E YVYY ١‏ 
متوسط ۲۳۹/۹۸ رت 1۷۸ 
زیاد YW/Y‏ ۳۳/۳ ۸۸ 
غیلی زیاد ۷۸۵ ۱۹/۸۸ ۱۰ 


نقشه طبقات خطر روش فازی سلسله مراتبی در شکل ٦‏ نشان داده شده است. ملاحظه می‌شود که به دلیل اثر 
بالای جاده و تراکم پوشش گیاهی. مناطق مرکزی. جنوبی و غربی پارک مناطق با خطر زیاد و خیلی زیاد 
آتش سوڑیٰ را به خود اخحتصاص داده و نیمه شمالی Sok‏ درجات کمتری از خطر را نشان oala‏ است. 
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شکل -٦‏ نقشه طبقات خطر روش فازی سلسله مراتبی در پارک ملی گلستان 


در روش رگرسیون لجستیک لایه‌های رستری موضوعی به‌عنوان متغیرهای مستقل و مناطق دچار حریق تاریخی 
به‌عنوان متغیر وابسته (با ارزش عددی یک) در مرحله واسنجی برای مدل‌سازی رگرسیون وارد نرم‌افزار Terrset‏ شد. 
معادله به دست آمده ب‌صورت معادله ٩‏ است. سپس پیکسل‌های خروجی به‌صورت طبقه‌بندی شده در ۵ طبقه 
استخراج که مساحت هر طبقه در جدول ۳ آورده شده است. ملاحظه می‌شود که تنها حدود ۱۱ درصد مساحت 
پارک در این روش پتانسیل زیاد و خیلی زیاد در بروز و گسترش آتش‌سوزی داشته و حدود ۸٩‏ درصد پارک خطر 


متوسط به پایینی در بروز و گسترش آتش‌سوزی دارند. 


سال دوازدهم ارزیابی پتانسیل خطر آتش‌سوزی با استفاده از .... ۱۸۵ 


۳,۵ = —14.6385 + 0.2D + 0.213527LST + 3.085010N + 0.0068215 — 0.020914TW + 
0.000279A + 0.002415H + 0.018821C — 0.002181TP (4) 


که مخفف‌های هر پارامتر مشابه معادله ۸ است. در معادله ۹ نیز با توجه به ضرایب اختصاص al‏ می توان مشابه 
روش فازی سلسله مراتبی تأثیر عوامل فاصله از جاده دمای سطحی و تراکم پوشش گیاهی را مشاهده نمود. با این 
تفاوت که به ترتیب عامل پوشش گیاهی QA)‏ دمای سطحی (۰/۲۱) و در رتبه بعدی فاصله از جاده (۰/۲) وزن 
بیشتری را کسب نموده‌اند. سایر عوامل نیز با ضرایب کمتر و دارای اثرات مثبت (مثل شیب ارتفاع جهت و انحنای 
دامنه) و منفی (عوامل TWI‏ و (TPI‏ در تهیه نقشه نهایی خطر مورداستفاده قرار می گیرند. 


جدول ۳- مساحت طبقات خطر آتش‌سوزی در روش رگرسیون لحستیک 


درجه خطر آتش‌سوزی مساحت (کیلومترمریع) درصد درصد تجمعی 
خیلی کم ۸ 0 0 
کم ۲۷/۸ VY/Y TTY‏ 
متوسط AAA ۷۵ MAY‏ 
زیاد ۷۰۰۳ AM‏ ۹۷۷ 
خیلی 35 ۱۹/۳۱ ۳/۳۵ m‏ 


متعادل‌تر در قسمت مرکزی پارک» مناطق با خطر کہ وسعت بیشتری دارند و مناطقی که در حاشیه جنگل‌های 
متراکم‌تر قرار داشته و cleo‏ سطحی بیشتر و رطوبت کمتری را تجربه نموده و در حاشیه جاده نیز قرار دارند. خطر 


پروز و گسترش بیشتری از حریق را OUS‏ می‌دهند. 
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شکل ۷- نقشه طبقات خطر روش رگرسیون لجستیک در پارک ملی گلستان 


۱۳۱ جغرافیا و مخاطرات محبطی شمارۂ سوم 


در شکل ۸ منحنی مشخصه عملکرد و سطح زیر نمودار دو روش مورداستفاده در ale yo‏ صحت سنجی OLAS‏ 
داده شده است. در روش es‏ سلسله مراتبی فازی» سنجه سطح زیر نمودار عدد /AL‏ و در روش رگرسیون 
لجستیک عدد ۰/۸۳ ملاحظه می‌شود که صحت پیش‌بینی را مشابه LIS‏ و سطح m de‏ توصیف می کنند. 
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شکل۸- منحنی مشخصه عملکرد در روش سلسله مراتبی فازی (راست) و روش رگرسیون لحستیک (چپ) 


ازانجاکه مقادیر سنجه مزبور مشه یکدیگر اسست. به مقایسه فراوانی پیکسل‌های مناطق حریق مرحله 
صحت‌سنجی در طبقه‌های مختلف در دو روش نیز پرداخته شده است که نتایج آن در شکل ٩‏ ملاحظه می‌شود. 
انطو کا شیرہ peed‏ لاق سام dogs AV. ae alla‏ یکبس هاش Bits‏ کاران خرن کار کی 
مرحله صحت‌سنجی را در طبقه‌های با خطر زیاد و خیلی زیاد پوشش داده است؛ درصورتی که این مقدار در روش 
رگرسیون لجستیک حدود ۲۲ درصد است. از طرفی روش رگرسیون لجستیک حدود ۶۱ درصد منطقه را در محدوده 


سال دوازدهم ارزیابی پتانسیل خطر آتش‌سوزی با استفادہ از .... AV‏ 


طبقه با خطر کم و خیلی کم توجیه نموده که با توجه به قابلیت اشتعال بالای منطقه از لحاظ پوشش گیاهی و دما 
منطقی 7 تست . 


7 روش فازی سلسله مراتبی ‏ روش رگرسیون لجستیک E‏ 


درصد فراوانی 


خیلی 3b; 3b;‏ متوسط کم خیلی کم 
طبقه خطر آتش سوزی 


شکل ۹- مقایسه درصد فراوانی کلاس‌های مختلف خطر آتش‌سوزی در دو روش مورد مقایسه 


-٤‏ نتیجه گیری 

مدلسازی و پهنه‌بندی خطر وقوع حریقء یکی از ملاحظات مخاطرات منابع طبیعی در هر منطقه است. در این 
مقاله به مقایسه دو مدل شامل رگرسیون لجستیک و مدل ترکیبی فازی سلسله مراتبی در پهنه‌بندی طبقات مختلف 
آتش‌سوزی و تعیین عوامل مؤثر در این مخاطره محیط زیستی پرداخته شد. جهت دسترسی به سوابق تاریخی 
آتش‌سوزی‌های منطقه از سامانه پردازش ابری گوگل ارث انجین استفاده شد. بررسی نتایج ضرایب عوامل ٩‏ گانه 
مورد بررسی در هر دو مدل OLS‏ داد نزدیکی به جاده دمای سطحی و ۱9۷1 از مهم‌ترین عوامل در بروز و گسترش 
حریق‌ها بوده است. یاتیش و همکاران (۲۰۱۹) نیز در تحقیق خود به نقش مؤثر NDVI‏ و دمای سطحی در بروز و 
گسترش حریق اشاره نموده‌اند. همچنین در مورد اھمیت عامل نزدیکی به جاده و عدم تأثیر بالای شاحص TWI‏ در 
نتایج مدلء پیائو و همکاران (۲۰۲۲) نیز به نتیجه مشابهی رسیده‌اند. در مورد دو عامل ارتفاع و شیب علیرغم وجود 
ضریب بالا در معادلات نهایی» ملاحظه شد که میانگین ارتفاع و شیب‌هایی که حریق در آن‌ها رخ دادم کمی بیشتر از 
میانگین منطقه است. این موضوع را منابع متعدد ازجمله پیائو و همکاران QI YY)‏ اسکندری و میسل (۲۰۱۷) و 
حیدری و همکاران (En)‏ و جانباز قبادی (۱۳۹۸) نیز عنوان نموده بودند. علیرغم اينکه هر دو مدل منحنی مشخصه 
عملکرد و سطح زیر نمودار مشابهی داشتند. بررسی پیکسل‌های مناطق دچار حریق در مرحله صحت‌سنجی OUS‏ داد 
که نتایج روش تحلیل سلسله مراتبی فازی منطقی‌تر از روش رگرسیون لجستیک است. چون حدود AV‏ درصد مناطق 


AA‏ جغرافیا و مخاطرات محبطی شمارۂ سوم 


قابل قبول عنوان نموده‌اند. ازآنجایی که در فازی نمودن لایه‌های مورد بررسی از تابع عضویت فازی خطی استفاده شده 
طبقات خطر در روش تحلیل سلسله مراتبی فازی OUS‏ می‌دهد که قسمت‌های مرکزی و جنوبی پارک که بیشتر جنبه 
تفرجگاهی داشته و دسترسی مناسبی از جاده‌ها دارند در طبقه‌های حطر زیاد و خیلی زیاد قرار گرفته‌اند. این موضوع 
لزوم فرهنگ‌سازی گردشگری» پایش مستمر مأموران محیطبانی و استقرار پایگاه‌های کمک‌رسانی در این مناطق در 


ماه‌های خرداد تا آبان را مضاعف می‌سازد. 
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